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Abstract 



We consider the zero-temperature dynamics for the infinite-range, non trans- 
\q • lation invariant one-dimensional spin model introduced by Marinari, Parisi and 

Ritort to generate glassy behaviour out of a deterministic interaction. It is shown 
that there can be a large number of met at ast able (i.e., one-flip stable) states with 
very small overlap with the ground state but very close in energy to it, and that 
their total number increases exponentially with the size of the system. 

1 Introduction 

A main issue in glassy systems is the analogy between glass-forming liquids and discon- 
tinuous spin-glasses, first pointed out in the pioneering works by Kirkpatrick, Thiru- 



malai and Wolynes ||KTW|| . In both cases the thermodynamical properties can be 



indeed related to the dynamical evolution in an energy landscape. In liquid theory one 
can define the notion of inherent structures ||SW|| (local minima of the potential energy, 
each one surronded by its attraction basin or valley) and configurational entropy, i.e. 
the logarithm of the number of these minima divided by the number of particles in the 
system. Then the low-temperature dynamical evolution can be described as a super- 
position of an intra-basin "fast" motion and a "slow" crossing of energy barriers. If the 
temperature of the system is small enough, namely less than the Mode Coupling criti- 
cal temperature T MC , the system gets trapped in one of the basins. Since the number 
of energy minima diverges exponentially with the size of the system, a thermodynamic 
transition can be associated with an entropy crisis: the Kauzmann temperature Tr of 
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the glassy transition corresponds to the vanishing of the configurational entropy. We 
refer the reader to [|MP|| for an overview on equilibrium thermodynamics of glasses. 



Consider now the class of discontinuous spin glasses, i.e. the mean-field models 
involving a random p-spin interaction. Also these models show a dynamical transition 
at a temperature T D (corresponding to T MC ) where dynamical ergodicity breaks down; 
a thermodynamic entropy-driven transition takes place at a lower temperature T\rsb 
(corresponding to Tr), at which replica symmetry breaks down with a "one step" 
pattern. Here the local minima of the free energy correspond to the solutions of the 
mean field TAP equations. Anyway at temperature T = metastable states with 
respect to any dynamics reduce to 1-spin-flip stable states. 

The main gap in the analogy between structural glasses and discontinuous spin- 
glasses is that in the latter models, unlike the former, the couplings between spins are 
quenched random variables. A significant, recent step in filling this gap has been made 
by the introduction of the deterministic, i.e. non-random, spin models which show a 
dynamical phase transition with a discontinuous order parameter and an equilibrium 
phase transition at a lower temperature associated with the vanishing of the high- 
temperature entropy ||MPK| , [BM] , |BDGU| , |JNM| , |PP|| . It is the high degree of frustation 



among the couplings, not the disorder, to generate a huge number of metastable states 
and thus the glassy behaviour. The discovery of these models proved that disorder is 
not necessary to reproduce a complex free energy landscape. 

Metastables states in infinite-range disordered spin-glasses have been extensively 
studied, both in the SK model where the number of 1-spin-flip stable states scales like 
exp(0.1992A) JrE], [BrM| , PGGQ , [MPVH , N being the size of the system (number of 



spins in the one-dimensional case), and in p-spin interaction spin-glasses |p^]| . 

Here we deal with the same question in deterministic models. By probabilistic 
arguments we will obtain, for the models introduced in |MP~RJ , a lower bound on the 



number of 1-spin-flip stable states, which increases exponentially with the size of the 
system. Hence this deterministic model exhibits the main feature of glassy behaviour. 
The paper is organized as follows: in Section 2 we review the basic properties of 



the model we consider, the sine model of ||MPR|1 , a deterministic, one dimensional 
chain of iV spins with long-range oscillating interaction. In Section 3 we study the 
limiting distribution of the rescaled energy density, showing that it gets (^-distributed 
in the thermodynamic limit. This property, which holds for the Curie- Weiss case, is an 
indication of the mean-field nature of the model. In Section 4 we deal with the state 
space of the system, computing explicitly (for iV prime analytically, for other values 
of N numerically) the distribution of the energy levels by flipping one spin at a time; 
among other things we show that there can be a large number of states with almost 
zero overlap with the ground state but very close in energy to it. Finally, in Section 
5 we derive the main result of the paper, that is a lower exponential bound for the 
number of metastable states at temperature T = 0. 

Acknowledgments: M. D.E. contributed to this paper during his visit to the 
School of Mathematics of the Georgia Institute of Technology, whose support and 
excellent working conditions he gratefully acknowledges. 
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this work was performed, for the kind hospitality during his visit under the exchange 
program between the Bologna and Brown Universities . 



2 Orthogonal interaction matrices and the sine model 

The basic setup is a probability space (Ejv, J 7 ^, P at). The sample space T, N is the 
configuration space, i.e. T, N = { — 1,1}^ whose elements are the sequences o = 
(ci, • • • , <7/v) with (jj = ±1; is the finite algebra with 2 2 elements, and the a 
priori (or infinite-temperature) probability measure P n is given by 

The Hamiltonian is the function on Sjy defined as 

H{a) = -\ J xy <r x <T y = ~< Jcr,cr > (2.2) 

xy 

where J = {J X y) is a symmetric real orthogonal N x N matrix given from the outset. 
Although many of the results presented here will hold for a generic symmetric orthog- 
onal matrix (e.g. of the form J = OLO T with L a diagonal matrix whose elements 
are ±1 and O a generic orthogonal matrix chosen at random w.r.t. the Haar measure 
on the orthogonal group) in what follows we shall examine the following particular 
example known as the sine model: 

J *V = , \ 8111 OAF , 1 • ( 2 - 3 ) 



V2N + 1 \2N+1 
which satisfies (we assume iV odd) [] 



N N 



JJ T = ld and X) Jr » = Z) J 2. = 1 - ( 2 ' 4 ) 



X=l 2=1 



This model has been introduced by Marinari, Parisi and Ritort as a deterministic sys- 
tem with high frustration (competiting interactions with different signs and strengths) 
able to reproduce the complex thermodynamical behaviour typical of spin glasses 
IMPllf . It has been investigated analytically in the high-temperature regime, (through 



an high-temperature expansion), and numerically also in the low-temperature phase 
(using Montecarlo annealing). The analytical study revealed the existence of a static 
phase transition at a temperature T$ = 0.065 where the high-temperature entropy 
vanishes, while evidence of the existence of a higher temperature Td = 0.134 where 
the system undergoes a dynamical transition of second order (i.e. with a jump in the 

^^Onc might also consider interaction matrices with zero diagonal terms, recovering orthogonality 
in large N limit. This amounts to put the average energy equal to zero (instead of —1/2) and may be 
convenient for particular purposes (see Section ||) . 
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specific heat) with a discontinuous order parameter has been put forward by numerical 
analysis. It has also been shown, using the replica formalism, that most of the ther- 
modynamical properties of this model are the same as those of a generic symmetric 
orthogonal matrix (the static transition corresponding to RSB while the dynamical 
transition being given by the so-called "marginality condition"). 



3 The limiting distribution of the rescaled energy 
levels. 

The knowledge of the eigenvalues of J imposes simple bounds on the energy of any 
spin configuration. Indeed a state vector a can be decomposed into his projections 
on the various orthogonal eigenspaces relatives to different eigenvalues. Here, due to 
orthogonality, the possible eigenvalues are +1,-1 so that 

- | < H{a) < |. (3.5) 

Let us consider the rescaled and shifted Hamiltonian (representing the energy per site, 
or energy density of the model, plus the 'zero point' energy 1/2) 

m = ^ + \ (3.6) 

which takes values in [0, 1]. We shall show that in the limit N — > oo the energy density 
h gets (^-distributed at x — 1/2. We point out that this property can be immediately 
proved for the Curie- Weiss model, thus indicating a mean field behaviour of the present 
model in the thermodynamic limit. To this end consider the partition function at 
inverse temperature (3: 

Z N (P) = e*p(-PH{a)) = 2 N E N (e^ H ), (3.7) 

where E N denotes the expectation wrt P N , and note that the characteristic function 
of h can be written as 

Bw(e -» ) = e -V»M^!0. (3.8) 

This expression will prove useful to compute the limiting expression of the characteristic 
function of the energy density h without knowing the expression of all its moments. 
To see this, we first decouple the spins as follows: let B be an orthogonal matrix 
such that B T JB = D with D = diag (d 1: . . . , d N ). Since det J ^ we have d, L ^ 0, 
% = 1, . . . , N, and det J -1 = Yl { d~ x . Let u G M N be such that a = Bu. We have 
< J a, a >=< Bu, JBu >=< u, Du >, and thus 

A N A 

exp ^2iv < ja,(r ^ = n exp (^v^ 2 ) 
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N 

n 



2,71 



exp 



-— + 
2 



_ ,_ U/'jOC'} I (XOC'i 



det J" 3 
(2vrA) iV / 2 



cxp 



< w, L> 1/2 a; > I cte 



cxp 



1 

2A 



< y, J l y > + < a, 



which, together with Q, ([Tj]) and yields 



Eiv(e 



-\h\ 



-A/2 



det J' 2 
(2yrA) Ar / 2 



exp(-^ < J x y > +^logcosh-^) (3.9) 



(As usual, the square roots appearing in the above formulas are only apparently ill 
defined: they disappear in the expansion because it contains only the even terms). 
The above integral can be evaluated by means of standard high-temperature expansion 
techniques which turn out to be considerably simpler if one assumes that ^ J a = 
(see ||PP|| ). As we have already noted, this assumption amounts to fix at zero the 
mean value of the energy. Also, the division by iV of the argument of the partition 
function leads to a convergence domain which is increasing as N itself. In this way, 
the asymptotic expression (for iV — > oo) of E n{g~ xh ) can be written in the form 



E N (e 



-\h\ 



-A/2 NG{\/N) 



(1 + OiN- 1 )) 



(3.10) 



where the function G{x) is an effective specific free energy. For the orthogonal inter- 
action matrix (12.31) one finds IPPII : 



G{x) 



Vl + 4x 2 - log 



1 + y/1 + 4x 2 



(3.11) 



It has the following expansion in the vicinity of x = 0: 



x 



G(x) = ^ + 0(x 3 ) 



(3.12) 



which, by the way, coincides with what one obtains for the SK model. This yields 



A_ 

AN 



+ o 



A^ 

AT2 



Summarizing, we have found that for any fixed A, 



Ejv(e- 



-\h\ 



-A/2 



1 + 0| ]v 



N 



oo. 



(3.13) 



(3.14) 



Using a well known theorem of probability theory [^J which says that a distribution 
function Gn converges weakly to G if and only if y27v(A) — > y?(A) for any A (where 
y?jv(A) and are the characteristic fcts of Gn and G respectively) and noting that 
y?(A) = e~ A / 2 is the characteristic function of the distribution fct G(x) = X[i,oo)> we 
then conclude that the distribution of h tends to X[±,oo)- 
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4 Flipping spins from the ground state and statis- 
tics of levels 



As already noted in [MFK], for special values of N the ground state, i.e. the configu- 
ration a G Sjv which minimizes the energy, can be explicitly constructed. Indeed, for 
N odd such that p = 2N + 1 is prime of the form 4m + 3, where m is an integer, let 
ex be the state given by the sequence of Legendre symbols, i.e. 



+1, if x = k 2 (modp), 
— 1, if x 7^ /c 2 (modp), 



(4.15) 



with k — 1, 2, . . . ,p — 1. Then (see the Appendix): 

H(a°) = 



(4.16) 



A typical ground state for p prime of the form 4m + 3 reflects the well known random 
distribution of the Legendre symbols (see Fig. [I], where a pair of ground states are 
shown for two different N values). No structure is present at any scale. Nevertheless, 
denoting by m° the specific magnetization of the ground state, i.e. 



m 



N N 

N - N 

x=l x=l 



X 



P 



(4.17) 



one observes that it tends to be a positive function of N, fluctuating around the 
value 1 / y/N. To let the reader better appreciate this fact we plot in Fig. ^] the 
total magnetization Nm° versus N. 



Figure 1: Ground state for N = 33 and N = 113 

For the remaining part of the paper we will restrict to p — 2N + 1 prime, with p = 
3 (mod4). We point out that this set has measure zero as a subset of the natural 
numbers. However, we have strong numerical evidence that some relevant properties 
that we are going to discuss hereafter, such as the behaviour of m°(N), the statistics 
of energy levels and the number of metastable states (see below) are somehow generic 
in N. 
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Let il s C Sjy be the subspace consisting of the (^) configurations obtained by starting 
from the ground state a described above and flipping exactly s different spins. Each 
point of Q s can thus be identified with a s-dimensional vector r e {1, . . . , N} s of the 
form r = (xi, . . . ,x s ), with Xi ^ Xj for i ^ j, which specifies the positions of the flipped 
spins along the chain of length N. We then define the 'flipping' map L T : Sjy — > S w 
as: 

In this way we can write 

n s = {L T a°} T . (4.19) 

The correspondence r — > a given by a = L T a° is plainly one-to-one. Therefore in the 
sequel we shall freely identify a state a = L T a° with the vector r. Alternatively, we 
can proceed as follows. Define the overlap q(a) of a given configuration a <E T, N with 
respect to the ground state cr° as: 

^) = ^X>^ ( 4 -20) 

so that g((T°) = 1. Then 

n s = {a E^ N : q(a) = 1-2^} (4.21) 

The following straightforward calculation yields the energy values on the space Q s : 
using the definition of L T , the symmetry of J and the fact that the ground state a is 
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an eigenvector of J to the eigenvalue 1 we have: 

H(Lro») = -^EE^X+^EE^X 

xer yer xer yg T 

+ \ E E J *»°l°l - 2 E E J -v a °A 

x£t y&r X^T ygr 

1 N 

= - 2 E J ^ <7 ° (j S + 2 E E J ^ <7 ° (7 S 
= -y + 2 E E _ 2 E E v^I 

x&t y=l x£t j/Gt 

= "Y + ^-2EEVW ( 4 - 22 ) 

xer yer 

Notice that for the Ising mean-field interaction: J xy — 1/N, one finds 

^, = 4 + 2.-2^=-^ (4.23) 

It is now possible to study the distribution of the energy levels on the individual 
subspaces Q s , where s = 0, 1, . . . , N. Let p s be the probability distribution restricted 
on Q s , i.e. 

Ps(C) =( N V E ( 4 - 24 ) 

^ ' aecnn s 

and let E s denote the expectation wrt p s . The n-th moment of the energy H on the 
subspace Q s is given by 

E s (H n )= [ H n (a)dp s (a)=( N ) Vr(l/), (4.25) 

so that the n-th moment E ^(H n ) of the energy on the whole configuration space S^r 
is 

e N (h») = jf FWdP.w 4e(!) E ^^ n )- ( 4 - 26 ) 

A tedious but straightforward calculation (see Appendix) yields the following expres- 
sions for the first two s-moments: 

N ( 2s x 2 



Vs(H)= -_M-_) , (4.27) 
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.; (g ).E,( g V(E.W)^ 4s(S -^:^ + Ar > (4.28) 



and consequently 

E N (H) = ~ E N (H") - (E N (H)Y = ^1. (4.29) 

These results indicate that, at variance with the ferromagnetic case where the energy is 
constant on each subspace Q s , here there is a significant overlap between the distribu- 
tions (for different s values) of the energy when restricted to fl s . In particular, from the 
espression of a 2 s we see that there can be a large number of states having small overlap 
with the ground state but nevertheless with energy very close to —N/2. For example 
we have cr%/ 2 — N/2, indicating that the energy restricted to the subspace may 
fluctuate over the whole energy range. This simple phenomenon is intimately related 
to the existence of metastable states and it will prove crucial in the understanding of 
the zero temperature dynamics, as discussed below. Fig. |3| shows the distributions of 
the energy restricted to various subspaces Q s . 




12.0 



Figure 3: Distributions of the energy over the subspaces f2 s . The size of the system is N = 23 
and s = 3, 4, . . . , 11 from left to right 

Another quantity of interest is the specific magnetization m(a) of an arbitrary state 
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a G Sat, given by: 

N 



N 

x=l 



In particular, given r = (xi, . . . , x s ), we have 

m(L r (7 ) = m(r) = m° - — V f - J . (4.31) 

Clearly, 



^ ] tGH s \i£t \PJ J ( ^ \ 25=1 ^ ' V ^ 



i.e. 



ren s 

s 

Moreover, we show in the Appendix that 



., > > 4s(iV — s) 4s(s — 1 



' x=l ^ ' 



iV 3 (iV- 1) iV 3 (iV- 

where d^(x) is the integer valued function giving the number of elements w G {1, . . . , N} 
such that u~ 1 x G {1, . . . , N} (where w _1 denotes the inverse mod p of u). As will be dis- 
cussed in the Appendix, d^{x) takes values around p/A with rather small fluctuations. 

Since YTx=x ff ) = ^' ^ ne ^ as ^ t erm i n ( |4-34 ) can beconsidered as a small correction to 



the constant value As(N - s)/N 3 (N - 1). 



5 Zero temperature dynamics and metastable states. 

We first introduce the following discrete 1— flip dynamics, given by: 

a(t + i) = / if H ^°) < H W> 

\ otherwhise, 

where, for each t, u(t) is chosen randomly in {1,...,N} with uniform distribution. 
Choosing an initial condition a(0) at random with respect to P at, one obtains a random 
orbit {cr(0), cr(l), . . . , cr(£)} for any realization {u(t)}i< t <e of length I. As a consequence 
of the previous analysis, we have the following remarks. 

• On one hand, it may happen that starting from cr(0) one reaches after t iterations 
a state a(t) G Q s , of the form a(t) = L T a° for some r = (x\, . . . ,x s ), such that 
H(L T a°) < H^L^LrCr ) for any to G {x±, . . . , x s }. 
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• On the other hand one can reach a(t) G Q s such that for some uj G {1, . . . , N}, 
L u a G n s+1 and H{L u a{t)) < H(o(t)). 

Due to the above observations, the overlap function q(a(t)) is in general not mono- 
tonically non- decreasing along a given random orbit (this at variance with the Ising 
mean field model). In particular there might be metastable states |[NS|| . Now, given 



uj G {1, . . . , N}, we shall say that a configuration o G S^r is uj -stable if 

H{L u a) > H(a). (5.35) 

Moreover, we say that a is 1-flip stable (or metastable) if it is cu-stable \/uj G {1, . . . , N}. 
We denote by n(N) the total number of such metastable states as a function of N. 
From ( 4.22| ) one readily obtains that 



H{L w a) = H{a) + 2 ^ J w a x a m (5.36) 
so that a is ^-stable if and only if 

{Jo) u (T u > Jw (5-37) 
Summing over uj and using ( [2.4D we see that if a is 1-flip stable then 

< Ja,a > > 1. (5.38) 



Recalling the expression (|2.2| ) of the Hamiltonian we see that a necessary condition for 
a G Sat to be 1-flip stable is that 

H{a)<V N {H) = - 1 -- (5.39) 

The main goal of this paper is to give an estimate of the number of metastable states 
for any given N. To this end we first performed some numerical investigations. For 
N < 30 we performed an exact enumeration of all configurations, whereas for larger N 
we run the zero temperature dynamics described above ( "deep quench" ) for a number 
or realizations {uj(t)} as large as 10 8 for bigger sizes, keeping track of the metastable 
states. As shown in Fig. ||, the growth of these states is exponential for generic values 
of the N. The best numerical fit yields 

n(N) ~ C ■ e°- 28N . (5.40) 

We remark that the same behaviour has been observed in ||PP|1 for the Random Orthog- 
onal Model; for spin glasses see |PP2 |. We now proceed to give a partial justification 



of this result by means of probabilistic arguments. 

Let t = (xx, . . . , x s ) and uj G {1, . . . , A^} be given. Using ( |5.36|) and Ja° = a , it is 
easy to see that if uj ^ r (i.e. L w a G O s +i) 



H(L w a) = H(a) + 2 ( 1 - - 2 £ J,.cr°a° 



1 
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Figure 4: The number of metastable states n(N) for N. The line represents the best fit 
n(N) f» e XN with A « 0.28 for values of iV such that p = 2iV + 1 is a prime of the form 4m + 3 
(filled points). Other points are for generic integer N. 



whereas, if uj e r (i.e. L^o" e fi s -i), we have 



H(L u a) = H(a) - 2 l + J uu -2^ J^>° ) 



If we define 



we then see that a configuration a = L T cr is w-stable if and only if 

h(r,u) < |(1 - Juu), ifw^r, 
h{r,uj) > \(\ + J w J), if w G r. 
12 



(5.41) 



(5.42) 



We now dwell upon the problem of characterising the behaviour of the function h(r, uj) 
so as the condition ( |5.42| ) can be effectively used to estimate the number of metasable 
states. Let us rewrite h(r, uj) in the form 

h(r,u) = ^= y %2&(u), (5.43) 



where 



U3X 



W) ■= [~ ) sin j . (5.44) 

Now, having fixed r and x G r, we can view the function £, x (u) defined in ( (5.44j ) as a 
random variable uniformly distributed on {1, . . . , N} and taking values in [—1, 1]. Its 
mean \i x and variance a x are easily computed: 



N 

1 ■ 



and, using (|2.4|), 



2 1 • 2 ( 2llUJx\ 



Here we want to study the behaviour of the sum r](r, uj) := Yl x &- £x{ u )- We remark that 
this sum, and thus h{r,uj) = 2r](r, uj)/ y/p, has to be regarded as a r.v. defined on the 
product of two probability spaces: for each fixed r, it is the sum of the i.i.d.r.v.'s £, x (uj) 
on the space {1, . . . , N} with uniform distribution (this comes from the very definition 
of the zero temperature dynamics); on the other hand, for each fixed uj, it can be 
regarded as a r.v. on Q s viewed as a probability space endowed with the distribution 
p s . Its mean is given by (recall that the symbol E s denotes the expectation wrt p s ): 

*.* - cr e z^rrt (:_-,>(.) (5,7, 

V ' t=(x 1 ,...,x s ) x£t V ' x=l V ' 

= -f(» = ^. 

N /_^^x\ , 2N 

x=l 

which does not depend on uj and equals s times \i x . Along the same lines one shows 
that 

Var ^ = w( 1 ~ji)- (5 - 48) 

Notice that unlike the means, here we have Var s ^ ^ s ■ a x . This discrepancy comes 
from the fact that, for any fixed uj, the sequence ^ X1 ,^ X2 , . . . , £ Xs is a sequence of distinct 
(and ordered) elements so that by no means we can view rj as a sum of independent 
and identically distributed objects. Nonetheless, also supported by strong numerical 
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evidence (see Fig. ^) , we claim that a version of the central limit theorem is applicable 
so that when N — > oo, s — > oo with s/N — > A, we have 

Ps ( a < < ^ - _L f e ~v 2 /2 dy (5 49) 



•v/Var s ?y / 

Assuming the validity of ( |5.49| ), performing the change of variables y = (x — X)/y/^y, 
with 7 = A(l — A), and setting a = A + a^/7, b = X+P^, we thus obtain an asymptotic 
gaussian distribution for h(- , uj): 

Ps (a < h(r,u;) <b)^ - 7 L= [ e-^-^l^dx. (5.50) 



Note that the r.h.s. does not depend on uj. One can actually say more: for any 
uj, uj G {1, . . . , N} we have h(uj~ 1 r, ujuj) = h(r, uj). Therefore the set of values of h( ■ , uj) 
on Q s does not depend on the choice of uj, i.e. {h(r, cj)} rgQs = {h{r' , u')} T , en , for all 

UJ, uj'. 




Figure 5: Distribution of the function h(r,uj) for a fixed u) G {1, . . . , N} and r varying 
in f2 s . Here N = 23 and s = 3 (left) and s = 11 (right). The solid line is the gaussian 
distribution in the r.h.s. of ( 5.50p . 



Having fixed an order for the lattice points (u\, . . . , ujn), ujj 7^ uj^, we now consider the 
following quantities: 



7r«,jv(wjfc) = p s ({> e fi s : H{Lu k o) > H{a)}) , 
the ^-probability that a randomly chosen state <r G Q s is u^-stable, 



(5.51) 



7r 8 ,jv(wfc + i|a;i,...,a; fc ) = p s ({a G O s : H(L Uk+1 a) > H{a) 

I H(L Wj a) > H(a), Uj = u t , . . .u k }) , (5.52) 
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the conditional ^-probability that a randomly chosen state a G Q s is c<jfc + i-stable given 
that it is Wj-stable for j = 1, . . . , k, and 

n s ,N = Ps ({er G fi s : H^L^a) > H(a), u = u x , ■ ■ ■ , uj n }) , (5.53) 

the p s -probability that a randomly chosen state a G £l s is 1-flip stable (i.e. stable for 



all possible flipping). Notice that by condition ( 5.42 ) the last quantity can be written 

as 

-1 s / i , T \ N-s 

I vv \ 



( N X e n o (V. *.) > n « (»<r, < . 

V / |. T , X»1=T 7 = 1 \ / 1=1 V / 

(5.54) 



(xi,...,X 8 )=T i=l V ^ J — 1 

(j/i,..-,yjv_ s )=T c 



The three quantities introduced above are related by the following identity: 

k s ,n = tt s (^i) • 7r s (uj 2 \uJi) • ^(wal^i, w 2 ) • • • n,{wN\wi, • • • , un-i) (5.55) 
and the total number of 1-flip stable states in T, N is, by definition, 

N /N\ 

We shall study the quantity n(N) in the thermodynamic limit: N — ► oo, s — > oo, with 

s/N — > A and < A < 1. In this regime we write 7r s jv = an d apply Stirling's 
formula to obtain 

with F(X) = -A log A - (1 - A) log (1 - A). (5.57) 



s 



y/27fN\{l - A) 



Note that F(X) is concave and symmetric around A = 1/2, with F(l/2) = log2. In 
this way we get for N large and s ~ XN with A ranging in the unit interval, 



N 



o 



2ttA(1 - A) 



log7T A 



n(N)~ / a — exp N [F(X) + ^— - dX. (5.5 



N 



It thus remains to estimate the probability tt\. Let us consider first the unconditioned 
probability Q5.51J) . According to ( |5.42|) and the total probability formula we have: 



7r s ,7v(u>fc) = ^-p s ( h(r,uj k ) > 1 



N \ 2 
+ (l - jj) Ps (h(r,u k ) < 



T 3 u k 



1 J(J k UJ k 



T^cok). (5.59) 



Here s/N and 1 — s/N are the probabilities that r 3 uj k and r 3 1 u) k , respectively. We 
can easily compute the conditional expectations 

fN- IV 1 2 . , s-1 fs-N\ 

E(/*|r9 W *)=( J _££&(u; fc ) = — U fcWfc , ( 5 ' 6 °) 
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and 

-l 



m\rjtu> k ) = ( N s X ) A^^ & ( a , fc ) = _^_(i_ J|iWb ). (5.61) 



In a similar way one can compute the variances 7 + (to>fc) and 7_(co>fc) conditioned to the 
events {r 9 c^} and {r ^ cc^}. For iV large and s ~ XN, retaining only terms 0(1), 
one gets 

E (h\r 3 u k ) ~ E (h\r $ u k ) ~ A, 7 _(o; fe ) ~ 7+^) ~ 7 = A(l - A). (5.62) 

Moreover in the thermodynamic limit specified above we write 7r S) jv(^fc) = ^x(^k) an d 
argue from ( |5.50| ) the following approximate expression for 

V 2?T7 y 1/2 V 27I "7 i-00 

1 + (I _ x] erf f (5.63) 



2 \2 J 
where we have denoted the error function by 



erf(z) = 4= / dx - ( 5 -64) 



It is not difficult to check that the r.h.s. of ( |5.63|) is convex and symmetrix around 
A = 1/2, where it reaches its minimum value equal to 1/2. 

Let us now come to 7r S) jv- In principle this quantity is to be computed by specifying the 
whole set of constraints embodied in Q5.54 ) or, which is the same, by computing the 
conditional probabilities appearing in eq. Q5.55p . However, this appears to be a difficult 
task. A first approach which drastically simplifies this task is to forget about the 
constraints implied by ( |5.54j ) and assume that (in the thermodynamic limit) the various 
cu-stability conditions become mutually independent, that is nx{wk+i\wi, ■ ■ ■ , uJk) = 
n\(uk+i), for all k = 1, . . . , N — 1, so that 



N 

tta = JJttaW- (5.65) 



fc=i 

Recalling eq. Q5.58Q and (|5.63| ), one is led to the following expression for n(N): 



"W"lV^A) eXP ( JVGl(A) ) dA ' <5 ' 66) 



where 

1 n .\ ./4-a 



G X (A) = F(X) + log \ + Q - A) erf (5.67) 
We show the shape of the the function G\(X) in Fig. || Evaluating the integral ( 5.66 ) 
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Figure 6: The function G\(X) of equation (|5.67| ) for A G [0, 1 



with the saddle-point method one gets 



n(N) ~Ce 



0.147V 



(5.68) 



Notice that the exponent is the half of what is observed numerically (cfr (|5.40| )). In 
the remaining part of the paper we shall argue that (|5.68| ) is indeed an estimate from 
below of the actual number of metastable states. 

The above discussion has been able to reproduce the exponential growth of the num- 
ber of metastable states with the size of the system. To understand the discrepacy 
between the estimated exponent and the one measured numerically one should note 
that the nature of the interaction makes the conditional probabilities play a major 
role in the asymptotics of the number of metastable states. To be more precise, our 
approximation which assumes mutually independent individual o;-stability events, i.e. 
71^(0^+1!^!, . . . ,Uk) — n\(uj k+1 ), is actually reasonable only for small value of k (this 
can be checked, for example, calculating the correlation functions). As numerical re- 
sults shows, for large values of k the specific form of the interactions make these events 
strongly dependent. In Fig. [7] we show the function P(k) providing the average of 
Tix{wk+i\^i, ■ ■ ■ ,^k) over a large sample of different permutations (ui, . . . of the 
lattice points. The conditional probabilities P(k) grow monotonically, almost linearly, 
from the initial (unconditioned) value up to a number close to 1. In other words, re- 
quiring that a large number k of spins produce an cu-stable state increases substantially 
the probability of doing the same for the remaining spins. 

Another way of understanding the constructive effect of the correlations is the following. 
Consider again the function h(r, u). Having fixed oJk+i-, we have already noticed that for 
s and iV large enough the values of h(r, u>k+i) with r G Q s are approximately distributed 
according to a gaussian probability density with mean A = s/N and variance 7 = 
A(l — A), regardless of the particular value of iOu+i- Thus, in particular, the same 
distribution are expected to arise if one considers the values of h(r, Uk+i) constrained 
to the subsets of configurations such that Uk+i G r, or uik+i t. On the other hand, 
if one picks u)\, o>2, . . . , 0Jk with Uj 7^ UJk+i, j = l,...,k, and computes numerically 
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Figure 7: P (k) versus k, for A/" = 15 (circle), 21(square), 23 (diamonds) and s = [N/2]. 
The conditional probabilities have been averaged over a sample of 10000 choices of 
(ui, . . .,u N ). 



the two conditional distributions of the values of h(r,Uk+i) given u^-stability, ... , 
a; fc -stablity (again with the constraints Uk+i G r or Uk+i 4- r \ one finds that their 
means move to opposite directions, thus increasing the probability of u; fe+1 -stability 
(see (15.59 )) . This is shown in Fig. where a system of size N = 21 and s = 10 
is considered. The two central distributions correspond to the unconditioned cases, 
namely the values of h(r, uj) for r G VL\q with the only constraints u> G r or u> r, 
respectively. Considering instead the values taken by h on the states r which, besides 
the constraints specified above, are stable with respect the first 10 spins, one finds two 
distributions whose mean values have moved towards opposite directions. An averaged 
over uj has been performed. 
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Figure 8: Graphs of the four distributions of the values of h described in the text. 



6 Conclusion 

We have investigated the statistical properties of energy levels and metastables states 
for a class of deterministic models, the most representative being the sine model [ |MPK| |, 
which have attracted much attention in recent years for their glassy behaviour despite 
the non-random nature of the interaction. We have pushed further on the analogy with 
glassy systems, proving a number of properties typical of disordered spin models. In 
particular, using number theoretic methods, we have described the energy (equivalently, 
free energy at T = 0) landscape as a function of configurations with a fixed overlap 
with the ground state. The analysis revealed the existence of states very different 
from the ground state but with energy arbitrarily close to it: this corresponds to the 
"chaoticity" property of spin-glasses systems, well established in long range models. 
More importantly, some of these states can be local energy minima (equivalently, 1-flip 
stable at T = 0). They are expected to have a significant weight on the partition 
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function in the low-temperature region, giving rise to the non-equibrium behaviour 
observed in annealing Montecarlo experiments. We have been able to estimate the 
approximate number of these energy minima. The analytic computations, combined 
with the numerical findings, strongly support the conclusion that the bound ( |5.68| ) 
estimates from below the number of metastable states n(N), proving their exponential 
increase with the size of the system. 

A number of basic questions about metastability arises now in a natural way, such 
as computing the energy density distributions of metastables states, studying energy 
barriers among them and their attraction basins. Stability of configurations with re- 
spect to the flip of an arbitrary number of spins is an interesting question as well. 
These problems are currently under investigation using the approach developed in this 
paper and will be addressed elsewhere. 



Appendix 



(Proof of ( 4.16 )) Choose N odd such that p = 2N + 1 is prime of the form 4m + 3, 
where m is an integer. Denote by cr° the spin configuration given by the sequence 
of Legendre symbols, i.e. 



+1, if x = k 2 (modp), 
— 1, if x 7^ A; 2 (modp), 



with k — 1, 2, . . . ,p — 1. Let us show that cr° = (cr° • • • o"5v) is the ground state for 
the sine model or, which is the same, that a is an eigenvector of J with eigenvalue 



1. For basic results of number theory used in the proof see, for example, ref. |[Ap| . 



2 

Vp _ 

2 1 / / i2nxy\ ( x 
y — exp — I — ex P 

Vpj^i 2l V \ p J \pJ \ p 

changing x i— > —x in the second summation 



i2irxy\ 



1 



N 



El i27rxy\ ( x 
exp ( ) 



x=l 



P 



P 



-1 



El i2nxy\ 



x=-N 



using multiplicativity of Legendre symbols: — 

' p 



x 



and the fact that ( — 

p 



— 1 if p = 3(mod4) 



-x 



P 



-1 

V 



1 

iy/p 



' N 

J^exp 

2 = 1 



i2irxy\ ( x 
p )\p 



-i 



Ei i2nxy\ 
exp ( J 



x=-N 



P 
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using the periodicity of Legendre' symbols : 



x + p 
P 



1 

i\/P 



N 



El i2irxy\ ( x 
exp 1 ) 



V J \P 



2N 



being = by definition 
1 / i2irxy\ ( x\ 

^\/p^ exp V P J \PJ 

using the separability for Gauss sums 



-( y -)t 

Vp \pJ t^i 



exp 



i2nx\ ( x 



iy/P \1>J ^ " V /' J \l> 
evaluating the Gauss sum 



( %2-nxy\ i x 

x ^ +1 exp V~r) \p 



' fV z)^P = < 



which is the desired property. 



(Proof of ([4.27|) and Q4.28| )) We sketch the basic steps of the calculation. Set 



a 



We then have 
For n = 1 we get 

e(ee^ 

T6fl s y£w 



a 



N 



+ 2s 



N-j 
s-j 



2 EEM =^ee 



x£t y£r 



? la"" 4 

u rGO s fc=0 



x£t y£r 



N N N 

x=l y^x,y=l x=l 

N N N 
£ £ C 2JxyCr x <J y + £(Cl - 02)/^ = C 2 (iV - 1) + Ci, 



whereas for n = 2 we have 

(\ 2 N N 

xdui y£w / 2=1 u^z;u=l 



N N 
,xj^z,u;x=l yj^z,u,x;y=l 
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N N N N 

+ ^2 °zJzy(?z<?y+ ^2 czJuyCuCy + c 3 J xz a x a z + ^ c 3 J xu a x a v 

yj^z,u;y=l y^z,u;y=l x^z,u;x=l x^z,u;x=l 



N 



N 



z=l 



+ C sJxx + C 2 JzuVzVu + C 2 J Z z + C 2 J, 

x^z,u;x=l 

N N N 

+ ^2 c 2Jxz&x&z+ C 2 J yz (Ty(T z + C 2 J xx + Ci J ; 

x^z;x=l y^z;y=l x^z;x=l 

= c 4 (N - 1)(N - 3) + 2c 3 (iV - 1) + 2c 2 (iV - 1) + c 1; 
which easily give the desired identities. 



N N 

E E ^ 

X7^z;a;=l y^z,x;y=l 



X yO ' x Oy 



(Proof of ( |4.34Q ) Let us first extend everything to the set {1, 2, ... ,p — 1} which, 
p being prime, is a number field. Here we can exploit the multiplicative structure 

of the field and of the 'extended ground state' cr° = ^|J, x = l,...,q, with 
q = p — 1. With slight abuse of notation we shall use the same symbols Q s , p s 
and E s to denote the corresponding extended quantities. It is immediate to see 
that m(a°) = and 



E s (m) 



« \ ^ q ^ 



0. 



In order to calculate the second moment we consider 

2 



7\E (E 



s J 

hS Cp(x) ( 



where, for any given r E Q s , r 2 is the collection, with multiplicity, of all possible 
products (all the operations are mod p). For example, if r = 

{x±, X2, x 3 } then r 2 = {x 2 , x\, x%, x±x 2 , x±x 3 , £2^1, £3X1, X2X3, £3X2}. Also, for any 
given z 6 {1, ... ,q}, 



c p (x) = ^2 {number of times x E t 2 } = jj{r \ u E r and w 1 x E t}. 
In particular, if = — 1 then -u 7^ m"" 1 ^, Vm = 1, . . . , q, therefore 

c p( x ) = J2 



u=l 



q-2 
s-2 



q-2 
s-2 
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If instead [ | ) =1, then there exists u such that u 2 — .r. i.e. 



c p (x) 



q-2 
s-2 



q-1 



(q-2) 



q-2 
s-2 



q-1 
s-1 



Putting everything together we get the following expression for the variance 
a 2 (m): 

2 



<r 2 (m) 



(I 



MM? 



2 [ s-2 



+ |(,-2,(^, +8 



q-1 
s-1 



q-1 
s - 1 



q-2 
s-2 



4s(g — s) 
q 2 (q - 1) 



We now turn back to our the original lattice {!,..., iV}. Again we can write 



N 
s 



tGQ 3 \x£t 



1 P 



N 
s 



x=l 



X 

C N {X) | - 

V 



In this case, however, the multiplicity function cn(x) can not be handled as 
easily as before. In particular, given x G {1, . . . N}, we denote by T(x) the set 
= {ui, . . . , Ud N ( x )} given by the w's in {1, . . . , iV} such that u~ l x G {1, . . . , iV}. 
The cardinality d^(x) of the set T(x) is a non trivial function of x. It is shown 

in Fig. |9| for 1 < x < N. If now = —1, then clearly 



c N (x) = d N (x) ■ 



N-2 
s-2 



On the other hand, if ^| j = 1 (i.e. -u 2 = x), then (note that either u G {1, . . . , iV} 
or -w G {1,...,JV}) 



c N (x) = (d N (x) - 1) 



N-2 
s-2 



N — 1 
s - 1 



We can then use these informations and write 
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Figure 9: The function cIn{x) versus x, for N = 551 and N = 933 




Finally, we have the following expression for the variance cr^(m) of the magneti- 
zation m over the space Vt s : 




from which one easily gets formula ([4.34 ). 
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